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Запропоновано аналітичні моделі гравітаційного поля та геологічного середовища для 
горизонтально-шаруватого середовища з кількома густинними межами поділу шарів у ниж- 
ньому півпросторі. Перша модель являє собою рівняння сили тяжіння, отримане з розгляду 
нормального потенціалу сили тяжіння у локальному околі поточної точки вимірювань у 
двовимірному зображенні. Для моделі середовища, на додаток до відомих у теорії інтерпре- 
тації потенціальних полів основних класів контактних поверхонь (Нумерова, Остромо- 
гильського і Страхова), введено новий клас контактних поверхонь - - клас Чорного. Наве- 
дено зразки теорем розділення полів для випадку кількох однозв'язних об'ємів і кількох 
шарів, що не перетинаються (два найпопулярніші початкові наближення середовища ) у 
досліджуваній моделі. Описано новий чисельний підхід до визначення початкового набли- 
ження густинного контакту і його асимптот у багатошаровому геологічному середовищі за 
допомогою кількох "контактних" алгоритмів. 
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Основное направление научного прогресса -- 
зто создание общих математических моделей 
в геологий, геофизике и в дальнейшем -- создание 
геолого-математических теорий зай 


Из докл. "Математическоег моделирование в 
петрофизике при комплексной интерпретациий 
даннькх ГИС, керна и опробования" 


ПАНГЕЯ 


Введение. В последнее время большое внимание уделяєтся вопросам переинтерпретации 
картографического материала, полученного по данньгм сьемок 30-летней давности, ввиду резко- 
го ограничения обьемов геологоразведки. Однако вследствие повьштенньх требований к точно- 
сти и зффективности интерпретации такая переобработка больших обьемов геофизических 
данньгх требуєт новьїх подходов к численному обоснованию методов, используемьтх при под- 
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готовке и обработке соответствующего картографического материала. Вполне логично обте- 
динить развитье математические конструкции с современньгюм интерактивньтм программньвюм 
обеспечениєем, дабь комплексно решить проблемьт оцифровьвівания, визуализации и истолко- 
вания геофизических данньх. Проблема истолкования данньгх, наиболее важная с точки зре- 
ния приложений (впрочем, тесно переплетенная на уровне первичной обработки с первьтми 
двумя), решаєтся в последнее время в русле перестройки аналитической и вьшщчислительной 
базьг теориий интерпретации потенциальньгх полей, обеспечивающей становление новой (тре- 
тьей) парадигмьгт в теории и практике интерпретации. 

Современнье методьт обработки геофизических данньх развиваются для создания цифровквіх 
аналитических моделей поля и геологической средь и вьісокоточньх алгоритмов, ориентиро- 
ванньх на применение в пакетах ГИС, особенно при отображений входной информации и 
результатов ее трансформации с возможностью сопоставления с априорной информациєей, так- 
же представленной в графическом виде (так назьваєемая визуальная интерпретация, рис. 1). 

Собственно, "послойная" обработка представляєт новьій зтап качественной интерпретации, 
в котором возможности представления информации в виде злектронного журнала измерений 
|Якимчик, 2010| можно использовать для предобработки информации в среде специализиро- 
ванньх пакетов общего и специального назначений. Примером такой предобработки могут 
бьть процедурьт вьщделения регионального фона, например, путем минимизации невязки меж- 
ду набором данньгх и некоторьїм полиномом в среде Майаб; предварительная фильтрация из- 
мерений в среде Ехсеї; нахождение особьгх точек поля методом обратной интерполяции поли- 
нома Лагранжа в Маріе; редукции поля с помощью алгоритмов В.Н. Страхова, реализованнь6"е 
в Руфоп, и т. п. Одним из наийболее удачньх примеров цифровой качественной интерпретации 
можно считать обработку значений поля силь тяжести в связке ПО Маріпбо - СогеїОтаму, пред; 
ложенньй в работе |Якимчик, 2010), что, отнюдь, не исключаєт их заменьт в обозримом буду- 
щем на ореп-5о0игсе аналоги. 

Коснемся, однако, иного направления в рамках детерминированного подхода к интерпрета- 
ции данньхх потенциальньгх полей -- разработки новьїх аналитических конструкций геофизи- 
ческого поля и геологической средьї. Относительно математических конструкций поля отсьмлаєм 
зайнтересованньтх читателей к работе |Якимчик, 2009|, а среди моделей средь авторам наийбо- 
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Рис. 1. Карта гравитационньх аномалий Украийнь в среде свободно распространяемой ГИС Сипапіит СІЗ. 
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лее близки плотностнье модели типа контакта, одну из которьх предлагаєм вниманию чита- 
телей. 

При разработке таковьїх конструкций ради повьшения точности интерпретации данньгх 
геофизики обязательно следуєт учитьвать условия адекватности постановок современньм тре- 
бованиям к получению исходньгх геофизических данньгкх и их обработке. В целом, суть зтих 
требований сводится к таким пунктам: 

е вьісокая тОЧНОСТЬ и технологичность алгоритмов и программ, их переносимость; 

е адаптация моделей поля к заданию исходньх измерений на коротких профилях; 

е адаптация моделей геологической средьт к ее сложному (нелинейному неоднородному) со- 
стоянию. 

Определение контакта в известньх моделях средьт Страхова, Нумерова, Маловичко |Чорний, 
Дубовенко, 2002| традиционно замьгкаєтся на вьічислений разности ундуляций некой теорети- 
ческой модели средьіщ от некоторого опорного ее злемента, от которого зти ундуляции не сли- 
шком сильно отклоняются. Подвох состоимт в том, что при зтом необходимо заранее знать не- 
которьге злементьтг геометрии (глубину, форму, ориентацию ) данного опорного злемента (мо- 
дельного тела либо контакта ), как начального приближения итераций. От вьрбора удачного 
начального приближения зависит, в целом, зффективность любого из существующих итера- 
ционнь5х алгоритмов |Старостенко, 1978|. 

Предлагаем вьіведеннье для горизонтально-слоистой средьт с несколькими плотностньтюми 
границами раздела слоев аналитические модели поля и геологической средьі, в которьігх зти 
опорнье злементьт рассчитьваются в процессе итераций. Модель поля представляєт собой 
уравнение сильг тяжести, полученное из рассмотрения нормального потенциала сильг тяжести 
в локальной окрестности текущей точки измерений в двумерном представлений. Размерь са- 
мой локальной окрестности существенно зависят от размеров самой тяготеющей области и 
точности аппроксимации самой нормальной составляющей потенциала |Дубовенко, 2011). 

Основное уравнение поля. Рассмотрим в трехмерном неограниченном евклидовом прост- 
ранстве Е 3) две области б; и б», вложеннье друг в друга так, что расстояние между их 
границами дб; и 96. удовлетворяет неравенству 5 (96 |,902)20 при б|0 6» єб». 
Рассмотрим в области 6 » начало координат и точки области 6 - б о б» С» або. 96 ги 


будем обозначать через 85 з (51, 52, 53 ), а точки замкнутой неограниченной области б "з 
1/2 


-Е)6, 20 "-д6, -- через Хе (Хі, Хо, Хз); |х -Е| є у (Х,8)-| У (хі - 6; )2 


Пусть область б. заполнена массами с плотностью 6» (5), а область 6 26,16» 
массами с плотностью 06; (5) (рис. 2). Обозначим через П (Х) единичную нормаль к поверхности 
96 ", внешнюю по отношению к 6 . Тогда нормальная составляющая силь тяжести равна 


и(х) 


пат аб ук) і 


где 5 (5) 2 6» (5) - 5, (5), 5є6»,5(5)- 0, 5є б о, / - гравитационная константа. Вообще 
товоря, силу тяжести описьваєт функция 


и (х)- | вгаді (х)| єї |/8(Є)| втайу г (х,Є) | Є, 2) 
| 


поскольку значения силь тяжести -- зто значения модуля градиента потенциала, а не со- 
ставляющие притяжения, в том числе и (1) |Дубовенко, 2009). 

Упростим вьтгражение (1), рассмотрев локальную окрестность 0 точки х при большом ее 
"Увеличений", Тогда с известной долей приближения |Дубовенко, 2011| участок дб « 0 пд6 " 
границь д6 " будем считать "неограниченной" плоскостью, т. е. дб з |Х: Х3 з Х3 (Хі, Хо) з 
- соп5і, (Х, Х,) Є 50 |, где 5 2 0 г 96 7". Расположив в зтой плоскости прямоугольную си- 
стему координат так, чтобьі ось 0Х3 бьтмла перпендикулярна к плоскости, получим 
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ат іівзіззчзвувззззезчу вена 


Рис. 2. Модель контактной средьг в классе Черного. 


п()є2, хе(хі,Хо,Хза), 52 (51,6256з), О'-опб, 
о «опо, дбаогодбт. 


Пусть в окрестности 0" функции 6 (5) - сопзі, 5 є б,со ,б, (б)єсопії, бе 0" пб,, 
и525г, - ба 6» -- собственное подмножество 0" . Фактические размерьт области О за- 
висят, прежде всего, от размеров изучаемой тяготеющей области С» и точности аппроксимации 
вьтражения (1) |Дубовенко, 2011|, т. е. в нашем случає 


ду" (хі, 51) ду" (хі, 1) 
дп(х) дхз 








 Ї 881) 48: - о | 


45 452 453 | «Є. (3) 


Если є » 0 достаточно мало, то О вьбрана удачно. Границу 96 » - 0 70 дб » вьібираєм доста- 
точно гладкой: 


96» ехо ху ежі (жу ук) (ху у жк ) є 25, 1, о1,2,3; 25) зда їжу є С (592 (5, ) |, 


При таких допущениях о среде получаем потенциал притяжения 


М (хг) Я. 9 дЕдЕ, дез 


Мазі дк 





где К з Їб. 
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Если массь, сосредоточеннье в 0 , генерируют аномалий и (Х| ), вьітянутье, положим, 
вдоль оси 0Х2, то их лучше всего измерять по профилям, параллельньм 05. Но формула (4) 
малопригодна для обработки таких измерений. Проинтегрируєем вьтражениєе (4) по 8», т. е. 
вьіщчислим интеграл 


где цевр-ху У 286, - Хо, Ма 83 о Хз, т ви? ми? 43 з дм. Делая замену переменной 














тіф У ,2 аа З 
интегрирования У - тієФ, фу - , Фагсіє -, 7" - 7 18 ФІ , получим 
со8? Ф і2 со5  Ф 
ю. тіф// сов? Ф ме | М 
(шум) ям | з Ї со8ФіФа---зіпФ з 

2 2 5/2 42 22 

Фр / сов Ф ФІ Фі 

2 1 
и ЗМ ЗМ 








и" ом ПРИ ОР Чо нн о ЕФ, ) 


то | (ц, М) - при 55 -9 со, 65 9 -со, и, следовательно, 





53 -Х3 
що охо, ек 45 48, 
П | (5з -хі) я (ба -ха) от 


где 5 -- проекция 6 » на плоскость 0Х| Хз, К - Їб. Положим Х3 з 0, 83 з 83 (5:|), а 8, «б, 
тогда 








(2) 
ьо 5 Ь 2 0) 
2 -хі)7з 
и (хо) -к| ав, | . в ока 1 ок 5 
а 40 (8р -хі) за а а (81-х) У (51) 
где о ш бю (рум і - б (51), ас86; 20 -- кривье, ограничивающие снизу и сверху 


тяготеющую двумерную область 5. 
Модели контактньх поверхностей. Рассмотрим несколько популярньх классов контакт- 
ньжх поверхностей, для которьх гарантирована единственность решения обратньтх задач. 


Геофизический журнал Мо 1, Т, 35, 2013 163 


Ю. И. ДУБОВЕНКО, О.А. ЧЕРНАЯ 


Класс А. Х. Остромогильского О5і м, (а, 0). Доказано |Остромогильский, 1970, что в 
зтом классе однозначно определяются переменная плотность 6 (5|, 53) и граница раздела 
сред Хз (5) по значениям внешнего поля. 

Граница дб раздела сред с переменньжми плотностями задается уравнением 


-П'кбха (51), 5Є (а, 6, 
- 6 
хз (5) о 2, є Га, 0, (6) 
причем 5 Х3з (а) - 0Хз (0) - 0, Хз (6; ) - - кусочно-непрерьвная функция и 
|8ха (х, / 2с(х, за) "9, ах сажд, а»б, 
|8ха (ху асо - ху 39, родах, «б, а»д. (9 


(1, о) (1, о) 
Подклассьг Обі (ю (2) 


и синклиналей соответственно. 

Класс В.Н. Страхова 5:09 (а, б) контактов рассмотрен неоднократно |Страхов, 1970; 
1974а, 19746; 1975а, 19756; 1976). Граница дб раздела сред с различньжми плотностями зада- 
ется уравнением 


класса 05:01: 0 (а, р) можно считать моделями антиклиналей 


-Пу жбХа (51 ), 6: єГа,б), 
Хз (5 уч (3) 
ПЗ», олені 
При зтом 0Х3(а)- 0, П|, 2 » 0, Х, (56:)20 -- кусочно-непрерьівна и удовлетворяєт на 


концах интервала Га, | условиям 


|5хз (5, )| за (8 за)", ахЕ, хажб, а»об, 


їни, 
|5хз (6, Д2б (6-61)77, б-бав кр, б»0, (9) 
принимая на Га, | конечное (счетное) число зкстремальньхх значений, так что множество 


ехіг2 (5) | непусто. Здесь 
га, БІ 


26 22 (с ) евирехіг2 (Є. ), -16 «Пі 7і -2 (с. ) -іпі ехіг2 (Є, ), -архй», (10) 
Га.6| Га.6| 


где с|,сь Є (а, |, аточки с|, с» могут бьть либо внутренними точками отрезка (а, 0| либо 
совпадать с его граничньми точками. В последнем случае при условий, что Й; «Й», 
сьоа,с»г-р. 

Если (а, р) -- бесконечньй промежуток, т. е. д 2 - оо, б 2 4 ее, уравнение границь (8) и 
условия (9) принимают вид 


2-2 (5)--Їи бо (5) 8Є(ОБо»бо) | 206)-би |єаз|ЕЇ"", 





О обі 


| 2(2)-й» | «62 | 277, В»о, (1) 


где П,Пг»»0,а таба)чне| ехіг вої, птоз ехіт 1 
І-бо 60! І-Єо 501 
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при сі, С» Є | -5о, 5о І, причем для определенности положим -25 «|, 2-2; 4 Пй,, О« П| «Па. 

Для контактной границьт раздела слоев с постоянной плотностью установлена единственность 
Г|Чередниченко, 1978 а; 1978 б). 

Пусть 6) -- горизонтальная полоса, ограниченная сверху и снизу прямьшми 2; з Пі, 
22й2»,0« Пр «Й, «же,ипусть 6 разделена на два слоя с плотностями 9; и 6 бесконечно 
гладкой кривой дб , звездной относительно бесконечно-удаленной точки, которая имеет 
асимптотой прямую 2 - й, Й| «Й « й», причем 


|а(2у-п|«К|Е|?, с»о, || збо, 2(5)єд6. (12 


Тогда, если значение й и величина скачка плотности Дф з 02 - 0; известнь, то граница 096 
раздела слоев по внешнему полю восстанавливаєтся однозначно. 
Доказано |Страхов, 1974 6), что решение обратной задачи о контактной поверхности класса 
(12) единственно, если известньт скачок плотности Дф и конечная дуга кривой 7 (5) из (12). 
Введем специальное обозначение для класса (12): СВ а, В (- се, со) , Класс (11) будем обозна- 
чать 51" В) ( -со, со), Очевидно, 


Ові1- 9 (а, бус8і 9 а, Б), СВОР) (оо оо) с 8) (ою, ов), 


54 (а, рус б» ОО (а є), (3) 


Проблема разделения полей. В практике гравиразведки наибольшее значение имеют 
задачи, в которьжх внешнее поле создают сложньєе распределения масс, например, п попарно 
непересекающиєеся обтемьі б ,К -1,2,.., п, или п попарно непересекающиєеся границьт раз- 
дела слоев различной плотности. В задачах со сложньвми распределениями масс единствен- 
ность решения обратной задачи верна лишь в довольно узких классах, т. е. при наличий боль- 
шой априорной информации об источниках. Теоремьт единственности обратной задачи гра- 
виразведки для сложньх распределений масс основаньт на теоремах разделения полей, которье 
позволяют свести решение обратной задачи по суммарному внешнему полю п обьектов (рудньх 
тел, границ раздела слоев) к решению обратной задачи для отдельньгх обьектов по полям от 
зтих обьектов. 

Пример 1. Пусть внешнее поле создаєтся массами с носителем в виде обьединения и 
конечньх попарно непересекающихся односвязньх обьемов б |, К - 1,2,...,П, с границами 
дб к , неимеющими попарно общих точек. Плотности масс в каждом из обьемов б , предпо- 
лагаются ограниченньми. 

Теорема разделения полей. Если известнь внешние относительно друг друга гладкие 
замкнутье конечнье поверхности 5 у, К -1,2,..., п, строго внутри каждой из которьх У 
находится одна и только одна компонента 6 | носителя масс, то по заданному суммарному 
внешнему полю |вне обьединения об»емов б у ) поля компонент б | определяются одно- 
значно. 

Из зтой теоремьі легко получить теоремьт для сложньх распределений масс, аналогичньт6е 
теоремам Новикова и Раппопорта -- Сретенского (для обьемов 0 у, звездньх относительно 
заданной точки |Чорна, 2003|, либо обладающих средними плотностями). 

Пример 2. Пустьв п-слоистой двумерной среде каждьгй из слоев б | имеет постоянную 
плотность 9. Границьтг раздела слоев дб у - - бесконечнье кривье, имеющиєе горизонтальньте 
асимптотьюш 7 з п и сей (пе »0, Пк 2 0) слева и справа (рис. 3). 

Теорема разделения полей. Если кажаая из границ раздела дб у содержит горизонталь- 
ньй отрезок дб у (5 - сопзі, ак 26 20 ), известнь проекция некоторой его части на ось Ох, 
ах охо, п прямне у; разделяющие границь слоев, то по заданному внешнему (при 2«0) 
суммарному полю Л9 с точностью 90 постоянньх определимь поля ДА |; отдельньх границ. 

Следствие |Гласко, Страхов, 1977). Если при зтом известнь положения асимптот 7 з Й і и 
2-й і и скачки плотностей Л9 і "Траницьт дк (К-1,2,..., п) раздела слоев восстанавливаются 
по заданному суммарному полю Л9 однозначно. 
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нозаочоачовеччочезоччччо чо зт я 


Рис. 3 К разделению полей от нескольких слоистьтх источников. 


Система непересекающихся контактньх границ и плотностей є (5) восстанавливаєтся един- 
ственньм образом |Филатов, 1974). 


Рассмотрим теперь решения нелинейного интегрального уравнения 





6 зоні 
зссюч |в ТОН цу ки 
зо Ме за) 
гле 2-7» (6), асхбхр -- "верхняя", 2-71 (6), асб хр -- "нижняя" границь плоской обла- 


сти 6 ,. Положим в уравненим (5) 9 (Х)-кК "МУ їх), 2» (6)-Н зсопзі, 21 (5)-2 (5), где 2 (2) 
принадлежит одному из множеств (например, (13), либо другому) допустимьтх функций (линеалу), 
тогда получим нелинейное интегральное уравнение по отношению к искомой функции 7 - 7 (5): 


б 2 2 
-ху7ЯнН 
"ТЕЗЕ 5 , з. 96 (14) 
у. "Кен "ра С) 





Его можно рассматривать как функционал на некотором линеале. Найдем его приращениє 
лу (Х,2)-9 (Х,2-п)-4 (Х, 2), где п (а)- п (б)-0 принадлежит заданному линеалу. В итоге 
имеем 


ь 
ло (х, 2) 2 Піп 


а 


(Е-хзн? 
(Е-х)2 ж (2)? 


6 зай 2 
фе іль арх р 


(Е-ХУ2 ж 22 








а 








27 - - 
щ і | о. СРУЧЕ | (5 а. Й Й М 
р С З з Шовех ва 


2 з12(Е)ч0 (5) (531 , 412 (8) (5) (51 
3| Е-ю? (БО ПОР 5-ю нн 





| "і п? (646. 
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Отсюда 
г 246) | (2-х)? ча2 (5)  ою2-22(5 
Ї з ПО 5е Па ; |З з Й (фен 
вес Ям 0 зо 507 І2 (59 НИ СОЇ а КОХ" за Кб 
у; 2 2 
2 (2)40 й 7 (2)--0(5)П -3 (6-х 
2 г ГА (5)--8 (5) П (8912 (5) (5) В (2) (5-х) Ї НЗ(БдЕ. а9 





зе (2-х? а (бу УРСР Р 


Алгоритмьі. На оснований (14)- (15) можно построить несколько способов восстановления 
контакта 7 - 7 (5), а 45 406, позаданному полю 0 (Х), сах «а, (а,р|с ІС, 41, теза, 0 | « тез |С, 4). 

1. Способ Ньютона. Считая известньми а,0,Н , разложим в ряд (14) относительно 2 (5), 
ограничиваясь первьгм членом. Учитьгвая, что 











2. (50Ю7зНО | «5-юЮбяг? | фоЮУНУ 10192 
91 | (5-ху9а22 0 (Е-ху? НО? ЦЕ-х 22 2 оби 
находим 
заданій 
-хугзнН 27 
9 (ж Пп (5-ю а Ї (9 д2(54Е 


зо (Е-ху? а? (5) 0 5 (Е-ху я (8) 


Отсюда, полагая АД 7 (5) 2 7 пут (5)- 2 п (5), имеем 











27 ) - Н 
Пра АНЯ МОМОТ ВР де дело 
а (57Х)7 ЖЙп" (5) з (5-Х)7 Ю2о" (5) 
0 аа) 
-| ; Ат (5342. (16) 


з (5-Х)? жа (2) 


2. Способ "контакт 1". Равенство (15) запишем с точностью до Й? (2), т. е. 





46 





27 - "б 
| І п 5у46- Пп х) 2 (5) ; 
З "ЕК НА (6) зо он З РАХ В ЛЬ СЕЇ 





р 5 2 
- - 0 п 
| (6-17 -І2 (2-8 (5) п (5 со. 
д 


(5-х)? Га (58 (5) (5112)? 
Сложим полученное тождество с (14), перенеся предварительно в (14) 4 (Х) в правую часть: 


р р 2 аа 
| 2205) повз Пп (6-Х) 27 (Є) феЕз 


з (5-Х)7 27 (8) з (ох)? Га (БУР СОЇ" 
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Ь 2 2 у; 2 2 
-Х)7 -|2 -0 п - Н 
«| (5 - | 2 (5)--8 (5) й | ЗЕ РИРУРУ Км (5 щі а 
з 1(5- Х)7 я 12 (2) 8 (5) п (2917) з. 5 Юна (8 


Тогда, обьединяя первь6й и последний интегральг, стоящие справа, получим 











27 - Н 
| о п (545- п С ху" я ; фе 
У т З а зо ес Є СЕРЕН ЇЇ 
р 2 2 
-Х)7 -Ї2 -0 п 
2 оз (5-Х0)7 -П2К88 (2) ПО 2 еуЧЕ. 


з (5-х? ж 12 (Був (5) (УР 


В зтом равенстве, как и раньше, - 14 8(Х)21. В итоге с точностью до Й і (5) находим 














р р 2 2 
27 -ху"зН 
Ї й п ч5а Па о у 95-909. (я 
з ЛЮ А 10) зо ЧЕК ФТ2СБНННСОЇ 
Полагая 2 ці (5)2 2 (5)ЯЙо (5), п20,1,2,..., м По (5)20, приходим к уравнению 
22 (5) -хугянН 
че - їн 452 Ла ООН водою. (8) 
злет нар 5) ак МОР УРА 0) 


Оно разрешимо при известном 706 (5):2| (5) 2 20 (5)-По (5). Так как Пи (5) 2 2рчо (5)- 2паг (5), То 











р р 2 2 
27п (5) -х)у7ЯнН 
о по ороаівве шось 1 ово 
ні- о бно а. (5-Х жі, (5) 
р 
22п (2) 
Ї ; гла (5) 45, а87 


з (5-х)? ад (5) 


откуда очевидно: алгоритмьт Ньютона и "контакт 1" формально совпадают только на первом 
шаге, так как 7; е 70, ана всех последующих шагах они формально различнь. 
3. Способ "контакт 2". Сложим величиньт (14) и (15) и с точностью до Пп" (5) получим 














і 22 (8) праве Со ЛИ 
з (2-х) жі" (5) з (5-ху? жа (ФужР ОВУ 
г Фою2-аб (5 г. (охузано 

- | 2 (бен п 45-90. 
з ПЕ-ху? а? (8 б. ог на 


Отсюда, обьединяя первьгй и последний интегральтг в правой части, а также полагая 


2 пз (5)2 2 (5)9По (5), п20,1,2,..., и По (56)2-0, ахехр, 
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приходим к уравнению 











р 
2 
Го о Пн (5) 52 
з (5-Х)? нап (5) 
-хю?зн? ; (8-х - 24 
- Па - ; 8 | - нрав 19) 
з. Кен бяга (5) 5-х)? ад (Б? 


Зто уравнение, как и (18), можно решать итеративно, если задано 70 (5): 21 (5) 2 20 (5) По (5). 


Уравнениєе (19) также можно дать в виде, аналогичном (187). Полагая П (5) 2 2 пу; (5)- 2 п (5), 
получаєм 


5 нії 2 
зна де Пи в цк) 
А о 





ь 
і 22 п (5) 
(8-х 22 (2) 





(5-х)? - яп (5) 
пн вена З 2 3 гад ча ФІ Б. 197 
з: ПЕОХУ ФаДЧ6Й 








ь 
І 22 п (5) 
з (2-х)? кад (5) 


А что, если определять не только контакт 2-7 (5), ачб«р,ной "асимптоту" Н з соп5і? 
4. Способ Ньютона. Примем во внимание, что 


9 ці83ЮТяні | 22 0 9 ,, бою?яні (0 он 
ді 


цеоху чо" (5-2 а3 дн (Е-хуза7  (б-х) ан 





и запишем 


р 
внанвнН а за ВНС 
(5-Х7ЯН 





ь РОДИ ь 
-ху7зН 22 (8) 
9 (а Пп С у з Й Ї ; 5 ; 
ж, «КЗ ВАС ІОЇ) 0 оз (ОХ ЗА (5) 
Отсюда, полагая 2 пу (5)7 п (5) 94725 (5), Нам ЯНо РАН, п20,1,2,..., м Д26 (5)20, АН є, 
асе«р, и задавая Н о й 20 (5), получаєм 
р р 


ь 

да 2Н ура 
Ї ем ДІ (2) 43 Ї 2 у АН а По 5 -. и 42-00. 00) 
з 57) ВАЛ КО) зе о На з Чех жаДСо 











Зто же уравнение в ой форме (непосредственно в "приближениях") примет ви 
друг Д А 

















р р 2 2 
22 (5) (59Ю79Н 
Го іа (5) «| Нана ее 
з (59-Х) аа (5) з (5-х) яН з  СБе Фа і5) 
22 (5) 2Н 
ДРУ 5 2 (5) 48 й у На дб, пе0,1,2,.., (207 
з (2-х)? ад (5) , (ох? На 
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5. Способ "контакт 1". В отличие от п. 2, здесь помимо определения 7, (5), 428), бу- 
дем находитьи Н,, 425 20. Для зтого рассмотрим функционал 


5 ана 2 
9 (х2,Нуа Пп (5 5, 
з (5-Х) 2 (5) 





и определим на линеале его приращение Л0 (Х,2,Н)-9(Х,24й;Н -т)-9 (Х,2,Н), где 
й (а)е п (5) 0. После несложньх преобразований в итоге получим 


р 
пера 





22 (5) 
) 


з (5-Х)? жа? (8) хдзн? 





| (вохуда(нат)? 00 (ох? (аз)? 
«| т т 
; (-х?ян? (Е-ху? жа? 


5 осо. й 5 пан 2 
і Ко З ню РР р -бс) 35-х)" - (Н -бт) зіва 
а 








Е-ю?зн і? 33 0 ІЗ(Б-ху? я (Н квт)? 


р 
(-х)2 о23 (Р) ації о РОЦЕНИР С ето цю 


з ПЕ? на (3 З Е-ю?н(анеу Р 





ПЗ (Є). | (21) 





а 


Сложим полученное тождество с точностью до й? (5) с равенством (14), перенесем 4 (Х) 
в правую часть полученного соотношения и сведем логарифмуеі: 

















р р 
27 Н 
1 наревн о Ной 
а КЕР ЗА ГО) АБУ БИ 
Ь ол 2 сад 2 6 ій 2 
«її б (5-7 (Н кт) і (5-Х)7 ж (24) Р (2-07 зов 
и (-хчн? (8-х? і? » (5-х) (2) 
откуда 
ь ь ь Й в 
27 он -Х"Я(НЯ 
І ес павуаені ее (5 М СЮ 986-909. 
, (2-х) ж 22 (8) , (Е-худ ян 5 (2-х? (а (Вчи (2) 


Полагая в зтом уравнений 
2 пні (5) 7 2 (5)9Йо (5), Нан Но ЯТо, П20,1,2,..., и Йо (5)Є0, то є0, ака хр, (22) 

и задавая Но, 206 (5), 425 20, получаем 
р 
фе й 


п 
з (5охуб ах, (8) 


р 


-| 22, (5) 
зла ЗНА) 


(2-х) зн? 


пі 








РОН, 
ін емів Ге 6-9 (Х). 03) 


раза со о 
-хузні"" 
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Так как Й ду (5) 2 2 пзо (5)- 2пч 5 Трої ЗНО дз 7Н ро: ТО уравнение (23) примет еще такой 
вид: 














р р р 
ен) (5-х)? 
Га ФЕН Но 469 й в а (х) 
(5-Х)7 вд (5) У Со З зо (5охуб жах, (5) 
р р 
27 ) 2Н 
Ї Що да (5)484Ї37оНан 4, П90,1,2,. 23) 


6. Способ "контакт 2". Сложим величиньі (14) и (21) с точностью до Й? (5) и получим 


























р р 
27 ха (НЯ 
Ге асо ЛОНА і ери С ою 
з НИ на 5) ЗЕ ЦО з | (Е-ху? ж (2яПу? 
7 (Е-ху-н2 1 (8-ху2-а2(9 1, 
Ї и 48 | з й (Ф4б, 
з КвеюЮ"-н | МІ 0 оо 9 од 
Отсюда, полагая (22) и залавая Ні Но, 2 27209 (5), а 25 хр, будем иметь 
Б р 2 2 
2 (5-Ю7ЯН 
прив ран са два Ла о авод00- 
(5-Х) ка (5) хуяНо з. 19780 асо (5) 
0 (Е-х)2-Н 0 (6-02 -22 (5) 
Ї тлавч| ОА? (фув. 24 


п 
, Пе-ху? ної з ГЕ? на (ФІ? 


Подставляя Й пл (5) 2 дао (5)- 2 пн (5), Тоні З Нора 7 Н дет получим то же уравнениє в "при- 
ближениях": 





Ь б Ь 
22 п (5) (Е-ху?зн? 

Ї ; 2 2 пч2 (5) 46 | а нтН Н по? 48-| іп п 

з (507кіп 05) ки Со о » (2-х)? ча, (2) 








р 


р 
22 п (5) 2Н , 
9(х Ї ; 3 да (5) 484 о пн 453 (247 





б 2 2 у; о зб 

(5Ю7УН к590 оп. (6) 
Ї зо (НвнсНо) 14 Гала (5)-2п (5317 8, 
з (5-Х"ЯНа) ІЕ- ху? чад (БУР 








паб0,1,2,... 
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Вьбор нулевого приближения. Пусть кривая 9 (Х), СС Ха имеет на ІС, 4 | несколько 
зкстремумов в точках с; є|с,41, і -1,2,...,П. Так как 4.» 0, то каждому зкстремуму отвечаєт 
одноименньй зкстремум 7 - 7 (5), а «5 «р. Определим во множестве і 900) найбольший 


и найменьший зкстремумь (наийбольший максимум 72с и найменьший минимум /|), т. е. 


95 -9 (сі )-5ир( ехіг 9(Х), 9; 29(с; )-іпв( ехіг 9 (Х)). 
с, 4) ге, 41 


Обозначим через 7с - 7 (Сі ) - 8ар | ЄХіГ 2 (5)), 2; 22 (С ; зві ехіг 7 (5)) соответствующие 
гс,4) гс,4) 

зкстремумьі 2. Здесь абсциссью 0 с и 76 (и 6; и 7;) совпадают, что не всегда имеет место, но 

в данном случаєе верно. Нулевое приближение контакта можно задавать, по меньшей мере, 

двумя способами (Ох направлена вверх): 


20 220(5)2--Но 44 (5,Но), ах ех, (25) 
20 220(5)2-Но -Код9 (5), ахбхр, (26) 
где Но 2»0 -- задаваемая константа -  "асимптота" функции 7 налевом конце промежутка 


Га,р1; 49 (5)- 9 (4)-9 (а); 49 (2, Н о) -- аналитически продолженное значение Д9 (5) на 
уровне -Н о; Ко -К(Н о) -- подлежащий определению козффициент. 

Рассмотрим, прежде всего, способьт определения Ко. 

1. Пусть в точке сє|а,р| располагаєтся 726 или 2;. Очевидно, 


; - оА2 2 Б 2 2 
дссуа п 5-01 ян 45- Га 091 -и, 
я кроку а 0) зі іезо Ко 





где йа85-с, ди-45, а, за-с, б, 2б-с. Возьмем 7с 22 (5), (2; 42 (8), ас5 хр, Тогда 





іо (5-3? я2? (Буре Е-с)7 нас Ь Ш(5-с)7 2 (ВУ Е-С)7 -20 1. | (27) 
Отсюда следует 
б 2 З 01 р. 2 
ЯН "ЯН 
де го, де ги. (28) 
ар 1 Ж25 ар М ЖЯі 
і: и?зНн? 
Оценим интеграл Пл --77--Фи, полагая, что а-0, рої, аро-с, б) 21-С. После 
шезьг? 
а1 
ряда преобразований окончательно находим 
бек 2 ц2 ЛЕТ 
ЯН 1-с)" ЯН УН 
Па" З" -фича тн сш? а 
оз а-су? жі? сіні? 
Н 
42Н агсів ц--- оо 27 агоїв ц. 
п с о 22 -5(1-К) 
Позтому 
а-с)чано сан? 
фено іно ЇЙ но 
азс)? жі сізає 
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2 
12Н агсів - Щ -215 агсів 9, (29) 
Н-с (1-0) 25 З Сі МЕЗО 
а-сруян? срян? 
др «ос всю з 
а-су)? жа? сужії 
Н 2 
ж2Н агсів Й -21; агсіє ; ; (30) 
Н 7-0; (7-С |) гі р (1-64) 


Полученнь6е неравенства являются исходньми для определения 7с и 7; при заданньх Н, 
д9с и 9). Отьскание значений 7с и 7; однотипно, позтому в дальнейшем индексьг опускаются. 
Составим функцию 


(а)за-да- су? ча? унсіп(с? кг? уч 2гагів оц 
22 -с(1-С) 


-а-сушіа-с?-Н21-сиооН2)-2н РЕ НН (31) 
Нігерії 
Вьшщчислим первую производную 
бозеньсу 52- офі офі с 
поесучяя?  осб4о 27 С(1-0 





27 або) 2 
- . с- 
"кре о 





2 


ПА 
2 7н6Йоб) 


2 
22 22 7 22|27 4 с (1-С 
-(ос) С 4 Загсіє -- ФУ 


оку за" | рожа 209 СЛЇ С) (2 04 Є ва 








(32) 


Для вьічисления Ї "(72) производную Ї "(2) лучше представить в виде 





, 2 7 1- 
Р (2)е2(-0с) зале удапиів Здаю -- 
(1-с)2 22 с? 7 7 


227 2с (1-02 





а яз , (327 
сао СУ) 820 ї20 ес(1- С) 0 


поскольку агсів ло зр агсіє зб агсіє 5 : . Дифференцируя второй раз (327 по 7, получаєм 
2 2 77 -с (1-С) 





рана 44-01? я 3 (322 вс (1-С)) 
і-субьотує сбз27) Феї) во" 
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ба (922 «20 (1-С)41) (27 ес (1-0) 
(22 -са-одр 22)2 





(33) 


Приближения 7; какого-либо из значений 7с и 7; определяем по схеме Чебьшева |Бе- 
резин, Жидков, 1960): 
2 2 
Р(2)  Р"(2а) Р (п) 
Р С з і (34) 
(а) ба) 





Последовательньюе приближения корня 210) уравнения Її (2)-0 можно получить значи- 


тельно проще, чем сделано П. Л. Чебьшевьм. А именно, если корень 20) располагаєтся в 
достаточно узком интервале | а, В|, то, если 7; вьвбрано с того же интервала, с точностью до 
мальжх 3-го порядка относительно 7 -7пі 100) 2 2, ЄГо,В| имеем 0- сс29) уз Р(Сп)я 


Р"(2 
зба 202 -2, У а о 


Отсюда при Ї "(2 )є0 (зто может бьіть, ибо только для кратньх корней Ї (2) 2 Ї ((2)-0) 
получаєм 


200) РО) 0 Р "(п ) 


Пи Ми 2 
п "ноу 2) та. 





Учитьвая малость разности (7 2 ра 78 у» , на оснований предьщдущего соотношения положим 


2 2 2 

ї ї ї 

сабо За «Ро готкуда 210) -2, 2- Ю а С отануві й Здьо 
Й "(2п) Р (22) 22)? 


Ребро баб (2) 
МИСРРИРУГУРИМ 


сходится (3--4 итерации для получения 7 (0) с точностью до 10 Зб что согласуется с резуль- 
татами, полученньжми в работе |Березин, Жидков, 1960)). 

Замечание 1. Из уравнений (29) и (30) следуєт, что Ї (25 )20,а Її (2і)20.А прине очень 
мальх Н » 0 имеем строгие неравенства: 








, т. е. получим (34). Процесс (34) при удачном вьвборе 72; бьістро 


ї(25)»0, Р(2і)«0. (35) 


Позтому корни 78, 2ц соответствующих уравнений Ї(75,9с )-0, ї(72ц ,9; )20 дают 
приближения 72с и 2і снизу и сверху соответственно, что может сказаться на "качестве" ну- 
левого приближения 20 5 20 (5). 

Замечание 2. Согласно замечанию 1 7с и 2; определяются не очень хорошо при больших 
Н, особенно, если значения 7с и 7; достаточно большие, Позтому дополнительно можно оп- 
ределить, если 59 - 0 (0) - 9(а)є 0, кроме 7с и 2|,заданной левой Н , еще и правую "асимп- 
тоту" Н " кривой 2-2 ( 25). Зто делаєтся по описанной схеме, только в вьграженимй Ї (7) вме- 
сто 9с или 9; берется 59 для определения Н б 


Вьчислив так для каждого заданного значения Н |, і - 1,2,...,П,левой "асимптотьт" искомой 
кривой 2 22 (5), а 45 «р, определим значения параметров 20), й Н Зб ПА ПЛ 
а з2) ні, пду22б)-но, АНіеНУ-ні, 


59 -9(0)-9 (а) ' 90 ах -95 -9 (а), о -9; -9(а), 


и найдем 


Є 
' 1 ' 
СРО бі 
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а ау ма к АН | 
зр по | б9иіо || і ва! 


Рассмотрим схему отьскания нулевого приближения для контактов, описьтваємьвтх функциями 


где У - 





из классов "антиклиналей" 05Ї ен (а,б) и"синклиналей" 05ї о (а, 6). Типичньми предста- 
вителями зтих классов являются функции (рис. 4): 
а) антиклиналь 


2(5)--Н за(1-86)8, 0«Ех1, (37) 


б) синклиналь 


2(вувен за (Пебув, КЕ СІ, (38) 


Для отьіскания Ко систему координат удобно вьібрать так, чтобьі точка зкстремума 9 (Х) 


находилась в начале координат, т. е. вано 0029 (0), тогда и интервал (|С,4| вьбираєтся 
"симметрично" относительно начала, так же как и а--с, б «С. При зтом для модели (37) 
(слева на рис. 4) имеют место соотношения а засіє С/Н,На| «242, 4, 2 02, адля моде- 


ли (38) -- соотношения й, загієС/72,На| » 202», б «2. 


Функция Р(2у«2сп (с2 з 22 уз Аг ага З - дсп (2 ЯН 2-4Н аб где 9-05 -- 
7 


в случае антиклинали, 9 - ді - - в случае синклинали (оба значения -- со своими знаками). Ее 
, С и 4с 
производнье Ї "(2)- 4 агсіє-; Р"(72)с- зі Так как Г(Н)- 9, то в случае синклинали 
2 с 37 
ї 4с 
(Н)«0, авслучае антиклинали Ї(Н )» 0.В обоийх случаях Ї"(Н ) зл зрету 0. Позтому 
сен 


РЕНУР"(Н)»0 длясинклиналимй Г(Н)Ї"(Н )«0 для антиклинали. В обоих случаях для оп- 
ределения 7с и 2; используется процесс (34). Попьтки обойтись грубьми приближениями 
25 й 2і, вьічиСсленньми не по (34), вообще говоря, успеха не имели!, хотя в ряде случаев 


приводили к хорошим результатам. ай 2 ; й 
Рассмотрим некоторье из попьток. Пусть вначале 0 «2СсІіп - --- - -4Н агсіє ні 4 7агсіє -. 
2 


ЧР 
Для синклинали (см. рис. 4) 2-Н чи, Н, 2,и»0, є, 2 (2, (знак равенства достигается 





С С 
при с -з о) и, следовательно, 4Н що ти дано заНнар-4(Н чц)а» 20, Позтому 
и 





Рис. 4. Модели контактов, описьваємьжх функциями "антиклиналей" и "синклиналей". 


| Слабая сходимость (до 50 итераций) и низкая точность, численная неустойчивость. 
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заці. З 29 -5 
іно ооо 9 3079 Но)є 20, АГ ЧО Фу 9-й. 0130) 
сао? Є 


Для антиклинали 2-Н -иИ,Н,и»0, а; 2 02 и, следовательно, 
С С 
4Н тині -) агсів- з4Н а - 420. 20. 
2 


Позтому 





9 

2 3 о 

УН 

і з бань бібванісо 2-2, (40) 
соні 96 


При небольших значениях и (по сравнению с Н) можно поступать по другому, а именно принять 





2Н 

шо" єн шу Телько 3 б) є 3-Й, 

с?зн? 

С С С СН 

(Н ки) агсіє -Н агсів - є | агсія -- -------- | Й 
НИ Н НО с2ан?2 
С 
для синклинали. Позтому 9 -- 4 окт и и, значит, 
С 
и -| 9 Ічаавко (41) 


поскольку для синклинали 9 « 0. Для антиклинали, хотя 9»0,но 25єН -й. Позтому формула 
(41) остается прежней. Заметим, что Г(Н)-9, Р "(Н )-4агсіє "пі и позтому (41) является пер- 


вьш приближением в итерациях (34). 
Второй способ основан на определений различньїх моментов кривой 7 з 2 (5). В простейшем 
оо 


случаєе, по теореме К. Гаусса Ї 9 СХ) дх «З2лт (т -- возмущаюшая масса, сосредоточенная в 
области, ограниченной кривьми 7; - 2(6), 22 - -Н ,ахехр,осьог-- вверх), позтому 


р р р 
те о 26) вч ондва | (2(8) нНУдЕ, оз. 
а а а 


Отсюда 75 з Рон, 7| «ТОН или 
р-а -а 
75 р и 90с04х-Н, 2і Р 9(х)дх-Н. (42) 
2п(-а) 2(р-а) 


со | 
В последних формулах полагаем Ї 9 (04х є Ї 9 (ху) дх. 





що с 
д д 
(Е-ху? ян? 
Обратим внимание на вьтчисление Ї 9 (х) дх з | іх Ї г - Я 45. Подставляя 
с бої З6О8Р РА (6) 
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цаф-х, ми2-бїх, цроб-с, й» 26-44, получаем 


4 р 2 аз Ш 
ЩО и и 














и з и" я7 о 
(вод ан (5-х) е-юЇ 
-| (2-х) -2Н агсів -22 (5) агсіє рака 
і (Е-Ху я" (5) Н 2(5) еаі 
й с?зн? д-еве С 
-| сі оВТи Ж2Н агсів -- 
й са (б) "РУ о з 


4 С а 
-2Н агсіїв - - - 22 (5) агсів - --- 4-22 (5) агсів - - -- дбс- 
сьо сі х| ; 





ооо с 


Із ід? «22 (6) -с1а Го? к22 (51-22 (Є) агоів -22(5) РН 


д 
2(5) 2(2) 


кача) сте? зН2у-Шасд?зн2)зон кову ЗОН ма | 


О численном решений. По предложенньм алгоритмам составленьт программьт и проведено 
численное моделирование. Теоретически зти конструкции показьвают несколько лучшую 
сходимость при вьчислений способом Ньютона по сравнению с классическими примерами 
|Гравиразведка ..., 1990|, но их практическая зффективность требует проверки на подготов- 
ленном материале полевьгх гравиметрических сьемок. Попьтки обойтись грубьшми приближе- 
ниями 7с и 7| , вьічисленньми не по (34), успеха не имели (крайне низкая сходимость бгіла 
на порядок меньше заявленной |Березин, Жидков, 1960)), хотя иногда приводили к матема- 
тически корректньтм результатам. Однако геологическое содержание грубьх приближений ед- 
ва ли достоверно. 
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